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Abstract 

The K~ optical potential in the nuclear medium is evaluated self consistently 
from a free-space K~N t matrix constructed within a coupled-channel chiral 
approach to the low-energy KN data. The chiral-model parameters are fitted 
to a select subset of the low-energy data plus the K~ atomic data throughout 
the periodic table. The resulting attractive K~ optical potentials are rela- 
tively 'shallow', with central depth of the real part about 55 MeV, for a fairly 
reasonable reproduction of the atomic data with x 2 /N ~ 2.2. Relatively 
'deep' attractive potentials of depth about 180 MeV, which result in other 
phenomenological approaches with x 2 /N ~ 1.5, are ruled out within chirally 
motivated models. Different physical data input is required to distinguish be- 
tween shallow and deep K~ optical potentials. The {K~ tov> ,ir) reaction could 
provide such a test, with exclusive rates differing by over a factor of three 
for the two classes of potentials. Finally, forward (K~,p) differential cross 
sections for the production of relatively narrow deeply bound K~ nuclear 
states are evaluated for deep K~ optical potentials, yielding values consider- 
ably lower than those estimated before. 
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I. INTRODUCTION 



There is a considerable interest in exploring the behavior of antikaons in nuclei and in 
dense nuclear matter [[!]] . The issues at stake concern the possibility of witnessing precursor 
phenomena, or even the onset, of kaon condensation in dense nuclear matter, as realized 
in heavy-ion collisions (see Ref. [[| and references cited therein) or in neutron stars (see 
Ref. || for a recent review). At present, the main evidence for a strong in-medium modi- 
fication of the KN interaction is due to the enhanced production of K~ mesons observed 
in subthreshold and near-threshold heavy-ion collisions in the KaoS experiments at GSI 
0-0. The extrapolation of the K nucleus interaction from standard nuclear density and 
zero temperature to the higher densities and temperatures which are relevant for the above 
phenomena is of course model dependent. Nonetheless, direct information on the K nu- 
cleus interaction at standard nuclear density and at zero temperature is extremely valuable 
and studies along these lines are therefore desirable. Unfortunately, even in this regime the 
situation is not clear cut at present. 

Systematic data bearing on the K nucleus interaction near threshold are almost exclu- 
sively limited to the strong-interaction shifts and widths of K~ atomic levels throughout the 
periodic table J7|. Although the 'atomic' K~ meson probes mostly regions of low nuclear 
density, the wealth of these data largely compensates in producing an effective constraint 
on the extrapolation to higher densities. The calculations existing todate for the K nucleus 
interaction at threshold essentially give two different predictions for the depth of the K~ nu- 
cleus potential at nuclear matter density. The phenomenological density dependent (DD) 
optical potential fits to the kaonic atom data HH, and the relativistic mean field (RMF) 
model calculations by Friedman et al. [[RJ, describe the data quantitatively very well, pro- 
ducing a deeply attractive potential (—Re K P t(Po) ~ 150 — 200 MeV). In contrast, chirally 
inspired models of the KN interaction, due to Weise and collaborators an d due to 



Oset and Ramos |L3]], give very good fits to the low energy scattering and reaction data in the 
strangeness S = — 1 meson-baryon coupled channel sector and to the K~p capture from rest 
branching ratios, but generally do not describe well the atomic data. Recently, following a 
suggestion by Lutz |Lj| , the K~ optical potential has been evaluated self consistently within 
such models [[L^,|l(J, yielding qualitatively reasonable fits to kaonic atoms |I7[] . Self consis- 



tency means that the outcome K~ optical potential should be accounted for in the in-medium 
K~ propagator that generates it within the appropriate scattering integral equation. These 
calculations predict a relatively shallow attractive potential (—Re V opt (po) ~ 40 — 60 MeV). 
Very recently, Baca et al. 1T8J improved significantly the fit to the atomic data by adding to 



the self consistent microscopic optical potential of Ramos and Oset |L5] a phenomenological 
s-wave 'tp' term of a moderate size (about 30% increase of the real part attraction, but 
50% decrease of the imaginary part absorption). However, this improvement was achieved 
at the cost of losing the direct connection of the optical potential to the chirally inspired 
microscopic model of the KN interaction. 

In the present paper, we aim at preserving the above mentioned connection by showing 
that reasonable parameters of the chirally motivated microscopic model of the KN interac- 
tion can be found such that the low-energy KN data plus the K~ atomic data are fitted 
simultaneously. This is accomplished in Section [HJ where the K~ optical potential provid- 
ing the atomic fit is constructed self consistently from the in-medium KN t matrix. This 
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in-medium quantity reduces to the free-space t matrix of the chirally motivated model as 
the density goes to zero. The calculations presented in Section |T| lead to a relatively shallow 
potential, Re V opt (p ) ~ —55 MeV. This depth is quite similar to that found in Ref. 



but the imaginary part derived in the present work differs substantially from theirs. The 
quality of the K~ atomic fit provided by our optical potential is superior. 

The vast difference in potential depths between the phenomenological and the more mi- 
croscopic approaches inevitably brings up the question whether at all, and to what extent, 
the K~ atom data contain unambiguous information on the K nuclear potential at nuclear- 
matter densities. Our experience leads us to believe that the behavior of the optical potential 
at densities about po is primarily determined by the functional form, or by the theoretical 
model used for the extrapolation from the low-density region which is more directly con- 
nected to the K~ atom data. Furthermore, a related question is whether there exists any 
experimental procedure for deciding which extrapolation is physically valid. This question is 
addressed in the present work by showing that the outcome of K~ initiated reactions at low 
energy is sensitive to the wavefunction of the K~ meson inside the nucleus, where different 
optical potentials produce noticeably different wavef unctions. As an example, we discuss in 
Section |T| the (fQ op ,7r) reaction into specific hypernuclear states, demonstrating that the 
production cross sections are quite sensitive to the K~ nucleus optical potential. 



Several authors have drawn attention to K~ nuclear states bound by a very deep 

potential, similar to the best-fit DD potential of Friedman et al. PH. The deepest states, 
bound by over 100 MeV, are blocked from decaying by the two-body mode KN — > 7rS 
and their width could then be reduced to about 10 MeV. If a deeply bound state of this 
kind is identified experimentally, it will provide evidence for a deep K~ optical potential. 
In Section [IV] we evaluate the (K~,p) forward cross sections for production of such deeply 
bound K~ nuclear states. Our evaluation yields considerably lower values than the estimates 



made by Kishimoto [19] for this reaction at 1 GeV/c. 



II. OPTICAL POTENTIALS FOR KAONIC ATOMS 

A. Empirical potentials based on KN amplitudes 

-^"-nucleus 'microscopic' optical potentials, i.e. potentials constructed from K~p and 
K~n scattering amplitudes as function of the local density, have recently been discussed 
by Baca, Garcia-Recio and Nieves in connection with deeply bound kaonic atom states 
fT8|. In this model the amplitudes are calculated self consistently within a chiral approach 



P~5|] and the optical potential is then constructed by multiplying these amplitudes by the 
corresponding proton and neutron densities. It was shown ]Tj| that such potentials yield 
marginally acceptable fits to kaonic atom data and that phenomenological modifications 
of the potential lead to much improved fits to the data. We begin the present section by 
adopting a similar approach. 

The interaction of a K~ meson with the nucleus in a kaonic atom is described by the 
Klein-Gordon equation of the form: 

V 2 - 2fi(B + V opt + V c ) + (V c + Bf] if, = (h = c=l) , (1) 
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where B is the complex binding energy and V c is the finite-size Coulomb interaction of the 
hadron with the nucleus, including vacuum-polarization terms. Equation ([!]) assumes that 
the optical potential V opt behaves as a Lorentz scalar. It is given by a '£(p)p' form: 

2pV opt (r) = -4tt(1 + ^)[a K - p (p)p p (r) + a K - n (p)p n (r)] (2) 

where M is the nucleon mass, p is the -^"-nucleus reduced mass, a K - p and a K - n are the K~p 
and K~n threshold scattering amplitudes evaluated at a nuclear matter density p = p p + p n , 
and p p (r) and p n (r) are the proton and neutron density distributions. The degree of success 
of any given potential in reproducing the experimental results for kaonic atoms is represented 
by the values of x 2 , defined in the usual way. However, comparing the values of x 2 achieved 
in the present work (see below) with those of Baca et al. one notes that our values 



are higher. The reason for that is most likely the fact that we use the experimental values 
for the yields of 'upper' levels || whereas Baca et al. seem to use the derived widths. It 
is easy to see that when the experimental accuracy is not high, the two definitions of the 
corresponding \ 2 may lead to very different x 2 values. The use of the measured quantities 
in the definition of the \ 2 function rather than derived quantities is the more appropriate 
way. 



Applying the amplitudes of Ramos and Oset |T5[ without any free parameter leads to 
X 2 = 300 or x 2 /N (x 2 P er point) of 4.62. This value should be compared with x 2 /N — 1-49 
obtained for a phenomenological potential || satisfying the low-density theorem constraint. 
The depths of the potentials at the center of a typical nucleus such as Ni are 44 and 54 MeV 
for the real and imaginary parts, respectively, compared to 187 and 71 MeV respectively for 
the phenomenological potential. In order to improve the fit to the data we followed Baca et 
al. Jl8| and added a phenomenological 'tp' term to the potential. Another possibility is to 
apply some scaling factors to the input amplitudes. In order to preserve the predictions of 
the chirally motivated model in the isospin zero sector where the strongly coupled channels 
produce the A(1405) resonance, we have applied scaling factors sr and si only to the isospin 
1 combination of the K~p and K~n amplitudes. More specifically, we have 



dl=0 — 2cbK-p — &K-ni a I=l — O-K-ni (3) 



and then we set 



a'i=i = $R Re a K - n + isi Im a K - n (4) 
and derive the modified scattering amplitudes as 

a K-n = a 'i=n a K- P = 2(^=0 + a 'i=i)- ( 5 ) 

Table | summarizes the results obtained using the input amplitudes, with or without 
adding a 'tp' term and with or without introducing scaling factors s. Also included in the 



table are similar results for the amplitudes of Schaffner-Bielich et al. \W\. Rather poor fits 
to the data are obtained when the potentials are constructed from the input amplitudes. 
The fits improve significanly when an empirical 'tp' potential is added, but the resulting 
total potential is not unique, namely, potentials giving similar values of x 2 may differ by 
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a factor of 2-3 in the nuclear interior. For the Ramos and Oset 
essentially confirm the depths originally found in Ref. 



fl5| input amplitudes, we 



181 upon applying this procedure. 



Further improvements in the fits are found when scaling factors, as defined by Eq. (|5]), are 
included. The values of x 2 /N may approach the best value obtained with phenomenological 
potentials pj but the required scaling factors are likely to be unphysical. 



B. KN scattering amplitudes in the nuclear medium 



Here we follow the chirally motivated model of Ref. |TT| for KN scattering and reactions 
near threshold. The A (1405) subthreshold resonance is generated in this model dynamically 
by solving coupled Lippmann-Schwinger equations for the tij elements of the t matrix in 
terms of the input chiral potentials v^. The six coupled meson-baryon channels included 
in this model are K~p, K°n, tt°A, 7r + S _ , 7r°S° and 7r~£ + . The use of a nonrelativistic 
formalism is justified for energies close to the KN threshold in accordance with the aims of 
the present paper. 

Applying the coupled-channel model to the calculation of the in-medium t(p) matrix, 
the density dependence of the K optical potential t^ N (p)p can be traced to the propagation 



of the A(1405) resonance in the nuclear medium. Waas et al. fl2"f , following Koch pT 



demonstrated that Pauli blocking of the intermediate nucleon is primarily responsible for 
the transition from a repulsive t^ N {p = 0), consistently with the subthreshold A(1405), to 
an attractive tx N (p) in the nuclear interior, as required by the A^-atomic data. Later, 



Lutz |14| stressed the importance of including the kaon self energy within a self consistent 
calculation, and Ramos and Oset |T5| provided a comprehensive calculation of this kind 
including also pion, nucleon and hyperon self energies. 

We extend the calculations of Refs. |f3~T| , [T~2|1 by including K and A" self energies. In this 
approach, the chirally motivated coupled-channel potentials are taken in a separable form 

Vij (k, k') = ^l^g^g^k' 2 ), 9j (k) = 1 , (6) 

Jtt 1 + l K / a j) 

where the momenta k and k' refer to the meson-baryon cm. system in the i and j channels, 
respectively, and the relativistic flux normalization factors $ are defined by 



with u>i, Mi and Ei denoting the meson energy, the baryon mass and energy in the cm. 
system of channel i. The coupling matrix CV, is determined by chiral SU(3) symmetry and 
includes terms up to second order in the meson cm. kinetic energies (see Ref. [[TT[] for more 
details). Finally, the parameter f n = 94.5 MeV represents the pseudoscalar meson decay 
constant, and the inverse range parameters were fitted to the low energy KN data in 



Ref. [111. Their values are 



OiK-p = ®K°n — 757.8 MeV , 

a n o A = 300 MeV , (8) 
a^+s- = avo^o = OV-E+ = 448.1 MeV . 
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The elementary amplitudes corresponding to the potentials (Bl) are of the form 



fij(k, k f ; E) = ^^ 2 )S#' 2 ) [(1 - C ■ G(E))- 1 ■ C sj , (9) 



ij 



where the meson-baryon propagator G(E) is diagonal in the channel indices i and j and is 
given by 

A ' P) ft E Ja i{p ) (2tt)3 kf -pi- H^, E hP] p) + iO ' 1 ] 

Here the integration domain Qi(p) is limited by the Pauli principle in the KN channels, p 
denotes the nuclear density and ki is the on-shell cm. momentum in channel i, such that 
uf = m| + kf, Ef = Mf + kf, and E = Ei + Ui is the total energy. In the denominator of the 
Green's function (|T0| ) we have included the meson plus baryon self energy Ilj. For simplicity, 
we neglect the self energy corrections in the pion-hyperon channels, expecting their effect 
on the KN channels to be secondary to the primary effect of including the K and nucleon 



self energies [16|. However, the pion self-energy effect was found nonnegligible in Ref. [15 
and this point deserves further study. 

The self energy term n in the cm. KN channels, close to threshold, is expressed in 
terms of the more familiar self energies of the antikaon and the nucleon as follows: 

n = ^H /? + ^n 7V , (ii) 

m K M 

where pkn is the KN reduced mass, = 2m^V^ t in terms of the K optical potential 
given in Eq. @, and = 2MV^ t in terms of the nucleon optical potential which in the 
present work was taken in the form 

K?t = , (12) 

where po = 0.17 fm~ 3 . For the nucleon optical potential we used Vb = (—60 — ilO) MeV. 
The real part is consistent with mean-field potentials used in nuclear structure calculations. 
The strength of the imaginary part is adopted from proton-nucleus scattering analysis |22 



It should be noted that when the meson and baryon self energies are turned off, our model 
reduces to the in-medium model presented in Ref. ||12|| . It is also obvious that the same 



applies for zero density. This means that we are able to reproduce the free-space threshold 



branching ratios and cross sections fitted in Ref. pl| . Below we show that a satisfactory 
description of these free-space low-energy data can also be achieved in a simultaneous fit of 
the model parameters to these data and to the i^~-atomic data. 



C. Fits to kaonic atoms and K p data 

The chiral model presented above was applied in x 2 fits to fT~-atomic data and to 
representative low energy K~p data. The latter consist of the three accurately measured 
threshold branching ratios [[£J 
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_ r(K-p^7r+S-) _ T(K-p -> charged) r(if-p -> 7r°A) 

7 r(ir-p -> 7T-S+) ' c r(lf-p -> all) ' n T(K-p -> neutral) ' 1 J 

plus four fCp - initiated cross sections at 110 MeV/c to the channels other than the 7r°A 
and 7T S° channels for which the quality of data is inferior. The adopted experimental values 



of the cross sections f25|-ftq] 



a(K-p -> K~p) , a(K~p -> K°n) , <r(if-p -> tt+IT) , <r(if _ p -> tt-£ + ) , (14) 

as well as the measured branching ratios 7, R c and R n are listed in Table [II]- The values 
calculated for the above seven quantitities are also listed, using (i) the original parameteri- 
zation of Ref. [0 (denoted 'no SC) which did not account for X and N self energies; and 
(ii) the present model (denoted £ Vn, full') which was fitted to both the free-space K~p data 
(0,0) as well as to the K~ atomic data, including self consistently K and N self energies. 
Table || demonstrates that extending the fit to include the atomic data does not spoil the 
good agreement with the experimental low energy K~p data. 

Table [TIT] summarizes the results obtained within the present approach. The upper two 
rows involve no fits to the atomic data; the K self energy is not included in the in-medium 
calculation pertaining to the first row ('no SC), whereas this self energy is included self 
consistently in the in-medium amplitudes of the calculation of the second row ('SC'). The 
nucleon self energy is excluded in both calculations. The \ 2 value for the seven K~p data 
points is of course the same in both calculations which, however, differ markedly in the Xatom 
values with respect to the 65 atomic data points. The calculated K~ optical potential is 
made considerably shallower and the description of the atomic data improves upon requiring 
self consistency, in agreement with Refs. [0,0 • The remaining four rows are for various fits 
to the 65 K~ atomic data points, some excluding ('atoms') and some including ('full') the 
7 free-space K~p data points in the fit, all with the K self energy included self consistently 
in the in- medium calculation. The last two rows ('Vn') are for the additional inclusion of 
the nucleon self-energy optical potential of Eq. ([121) . The table also specifies the resulting 
scaling factors which multiply the pseudoscalar meson decay constant f n = 94.5 MeV and 
the inverse range parameters a« of Eq. (§). It is seen that the modification of the free- 
space model parameters is moderate. Also given are the depths of the real and imaginary 
potentials at the center of the Ni nucleus. The depth of 55 MeV for the real part (last 
row in the table) is very close to the corresponding depth found in Ref. |T3, but the depth 



of 60 MeV for the imaginary part is considerably larger than that found there. It is clear 
from the table that, with relatively small modifications of the chiral- model parameters, it is 
possible to achieve reasonably good fits to both atomic and free-space data. The best Xatom 
means Xatom/^ °f 2.2, a value which is significantly lower than the value 2.7 corresponding 
to the hybrid fit of Ref. [0, yet considerably higher than the value 1.5 obtained with best- 
fit phenomenological potentials 0. It is therefore again concluded that when a free-space 
interaction model is used to fit the K~ atomic data, the quality of the fit is inferior to that 
achieved with phenomenological fits which are constrained only by the low density theorem. 

Figure [I] shows cross sections calculated for six free-space K~p initiated reactions in 
comparison with the data (as compiled in Fig. 1 of Ref. |l I|). The results of our full model 
('Vn, full') and the free-space calculation ('no SC) were chosen for illustration. The cross 



sections calculated for the other parameterizations of Table III are of a similar quality. It is 
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emphasized that we obtain equally good fits to the cross sections as well as to the threshold 
branching ratios even when only K~ atomic data are used in the fit. Moreover, we find 
that the inclusion of the nucleon self-energy optical potential Vj^ t plays only a marginal role 
in the improvement of the fit to the data. Of course, the size and form of V^ t influences 
considerably the depth of the resulting K~ nucleus optical potential, as demonstrated in the 
table. 

Figure |2] shows the density dependence of the isospin averaged (effective) threshold scat- 
tering amplitude a e s = (3a^g 1 +a^g°)/4 for three cases selected from Table |Hl|: (i) no medium 
effects beyond Pauli blocking are included ('no SC, dashed line); (ii) the self-consistent cal- 
culation including the K self energy ('SC, dot-dashed line); and (Hi) plus including the 
nucleon self energy ('Vn, full', solid line). The change of the sign of Re a e s from negative to 
positive corresponds to the transition from an apparently repulsive free-space interaction to 
an attractive one in the nuclear medium. In the 'no SC model, in which medium modifica- 
tions are represented only by the Pauli blocking effect, the transition occurs at p ~ O.lpo- 
When the K~ self energy is taken into account, this transition occurs at a lower density, 
p ~ 0.05^0) and the inclusion of V™ t pushes this transition density even further down. Nev- 
ertheless, the free-space (p = 0) threshold scattering amplitude remains negative even in 
this case, reflecting the dominance of the A(1405) 1 = subthreshold resonance. We note 
that since the 'no SC and 'SC parameterizations coincide, the low-density limit of a e g for 
these cases is the same, whereas a e g(p = 0) for the 'Vn, full' parameterization assumes a 
different value. 

Taking the K~ self energy into account generally leads to a weaker density dependence of 
the threshold scattering amplitude, both for its real and imaginary parts. This indicates that 
the K~ optical potential evaluated within such self consistent models is well approximated 
by a t cS p form (where t eS = const.) over a wide range of densities. A genuine p dependence 
of t e s appears only at very low densities. 

The free-space (p = 0) scattering amplitude //=o, as a function of the cm. energy E, 
is shown in Fig. [3[ The calculated amplitude for the 'no SC parameterization (first row 
in Table III, dashed line here) is compared with that for the 'Vn, full' parameterization 
(last row in Table [TTT], solid line here). The peak of Im fj =0 is shifted upward by over 10 
MeV when the K~ atom data are included in the fit. These amplitudes, for which the real 
part changes sign at the energy where the imaginary part peaks, provide a signature of 
the A(1405) subthreshold resonance. We note that a precise reproduction of the A(1405) 
spectral shape requires a more involved calculation than the 1 = KN — ttTj coupled channel 
calculation reported in the present work and, therefore, this spectral shape was not used 
here as a constraint. A new datum that was unavailable to the authors of Refs. [0,0 is the 
K~p scattering length, deduced from the recent measurement of the 2p — > Is X ray in kaonic 
hydrogen |2(J and which we too have not included in the fit to the data. Model 'Vn, full' 
does very well with respect to Re a K - p , but does poorly for the imaginary part. However, 
the rapid variation of Im /j =0 near threshold should make it fairly easy to reproduce it by 
slightly varying the parameters of the present chiral model. 



For completeness, we also tested the 10 channel chiral model of Refs. |[13jl5 |, including 



also the rjA, rjY, , K + E~ , and channels. We treated it on the same footing as our 6 

channel model, namely a separable form was used for the coupled-channel potentials, and 
only K and N self energies were considered. Fitting to the K~ atom data led to \ 2 values 
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comparable to those for the 6 channel model. 



III. STOPPED K~ REACTIONS AS A TEST OF THE K~ NUCLEUS OPTICAL 

POTENTIAL 

In this section we discuss the possibility of testing the K~ optical potential at threshold 
by studying (K~ top) 7r) reactions to specific A hypernuclear states. The (K~ top) n~) reaction 
has been used to explore the spectroscopy of A hypernuclei (see Ref. |?7| for a review). 
Recently, the use of the complementary (i^st op , 7r°) reaction was proposed p8] and first results 



on C are forthcoming [F2UJ. Calculations of A-hypernuclear formation rates for (K stop ,-K ) 
reactions were presented by several groups [pOj -33] within the framework of the distorted 



wave impulse approximation. The present calculation follows the approach described in 
Ref. |31| . Our aim is to study the sensitivity of the capture rates to the choice of the 
K~ nucleus optical potential, provided the latter was fitted to the K~ atomic data. Here 
we limit the discussion to the following K~ capture-at-rest reactions on 12 C: 

K- + 12 C — 7T- + \ 2 C , K- + 12 C — tt° + \ 2 B . (15) 

The capture rate per stopped K~ from the 12 C initial state i to a A-hypernuclear final 
state /, Rfi/K~, is given by |U| 

tf-.Bw^^ 1 " 11 '' , ,16) 

q f uj f 4np N 

where the fractions _R(7rA) are the in-medium branching ratios for K~N — ► 7rA capture at 
rest, as given in Table I of Ref. ||31|| , and the kinematic factor involving the pion momentum 
q and energy u> in front of R(nA) is due to considering the pion final-state phase space 
with respect to an average closure phase space marked by bars. The product of these two 
factors for capture on 12 C is 0.074 (for i?(7r~A)) x 1.456 (for qfUf/qftjjf), equaling 0.1077 
for (iQop, 7T~) and half of that, by charge independence, for (K^ op , tt°). The distorted wave 
(DW) transition amplitude is given by 

F £ W (9/) = / d'r X { q]\r) p fi (r) * nLM (r) , (17) 

with pfi denoting the nuclear to hypernuclear transition density matrix element. The brack- 
ets (• • •) in Eq. ([IB] ) stand for averaging on the initial substates and summing over the final 
ones. Finaly, pN denotes the effective nuclear density available to the capture process, 

Pn = J d 3 rp N (r) | y nLM (r) | 2 , (18) 

where Pn(t) is the appropriate nucleon (proton or neutron) density and ^ n LM( r ) is the 
K~ wavefunction in the nL orbit from which the K~ meson is captured. 

The radial wavefunctions for nucleons in 12 C, and for the A hyperon in \ 2 C and \ 2 B, 
were generated by solving the Schrodinger equation for a real Woods-Saxon potential with a 
diffusivity parameter a = 0.6 fm. The depth Vq was adjusted separately for each baryon orbit 
to yield the observed binding energies. The proton Coulomb potential due to the nuclear 
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core charge distribution was also included. The single-particle wavefunctions thus obtained 



were used in the calculation of Pfi(r) for the DW amplitude (|17D and also for constructing 
Pn(j) in Eq. flTEl). The value of the radius parameter r was chosen such that for protons the 
r.m.s. radius of p p was equal to the r.m.s. radius of the charge distribution after unfolding 
from the latter the proton size. As already observed in Ref. |[31||, since the baryonic radial 



wavefunctions that enter the capture rate calculation are nodeless and real, the resulting 
calculated rates flTE| ) are considerably less sensitive to variations in these wavefunctions than 
to similar model variations of the pion and K~ wavefunctions. 

To generate the pion DW Xq m Eq. ( |T7D we used the measured pion elastic scattering 

angular distribution for 6> cm < 90° at 162 MeV on 12 C []34| in order to fit the standard 
pion-nucleus optical potential 



2p n V: pt = 4tt[-(1 + ^)V(r) + (1 + ^)" 1 coVp(r) ■ V] 



(19) 



for which \ solves the Klein Gordon equation. We have managed to improve considerably 
the fits of Ref. pi , getting as low as x 2 /N = 3.2 for the following values of parameters: 



-0.24 + i 0.18) m; 1 , c = (0.20 + i 0.31) m^ 3 . 



(20) 



In addition to the chirally motivated K~ optical potential discussed in the previous 
section, several K~ initial-state wavefunctions were generated using the density-dependent 
form of the K~ optical potential |§ 



2pV c 



opt 



-Air 



1 + ^ 

M' 



a + B 




P{r) , 



(21) 



where the notation follows the 't(p)p' form, Eq. (0). For these optical potentials, the energy 
shifts and widths of the 2p and 3d levels in kaonic 12 C, and the corresponding wavefunction 
^nLM(^), were then obtained by solving the Klein-Gordon equation ([!]). Four different 
optical potentials, ordered according to their central depth, are listed in Table [TV]. The 
'chiral' potential corresponds to the relatively shallow potential of the present work. The 
deep potential 'DD' was obtained by Friedman et al. || fitting the parameters B and 7, with 
the value of a held fixed at the empirical K~ N scattering length, so that the potential (|21~|) 
satisfies the low-density limit. The potential L t e s ' was obtained for B = as the best-fit t e gp 
solution for the standard density independent version of V op t. Since the above potentials 
were obtained within global fits to the available i^~-atomic data, they do not necessarily 
reproduce precisely the experimental data on kaonic 12 C ||35|| , and for this reason we added 
in Table |TV| another potential ('t e fr') designed primarily to fit these latter data. 

The calculated capture rates per K~ , Eq. flTBl), are shown in Table |V]for the production 
of the 1~ hypernuclear ground states off 12 C, assuming atomic capture fractions f p = 0.23 
and fd = 0.77 according to a cascade calculation by Batty which fits the absolute and 
relative X-ray intensities for 12 C. It is clear that the deeper the K~ optical potential is, 
the lower the calculated rate becomes. This pattern is caused by the strong-interaction 
bound d state generated by all but the 'chiral' potentials, as may be recognized by the 
repulsive shift (see Table |IV|) they impose on the 3d atomic state. By orthogonality, the 
atomic wavefunction acquires then extra nodes within the nucleus, thus causing substantial 
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cancelations in the DW amplitude ([T7J). This effect was extensively studied in Ref. |3T] for 
K~ capture from atomic p states in 12 C, but for the dominant capture from atomic d states 
there was no similar effect since the relatively shallow potentials considered there (similar 
to the present 'chiral' potential) did not produce strong-interaction bound d states. All the 



calculated rates shown in the table are lower than the measured values. Tamura et al. [27 
report a rate of (0.98 ± 0.12) x 10~ 3 per stopped K~ for the production of \ 2 C. This value 
is still four times larger than our highest calculated value, that for the 'chiral' potential. 
The largest uncertainty in our calculation is due to the pion distortion effects which we 
estimate as about 10%, so the discrepancy between experiment and calculation cannot be 
resolved at present. Therefore, although the K~ capture from rest reaction does exhibit 
significant sensitivity to the type of V op t, it cannot yet be used to exclude potentials which 
are permissible from the point of view of fitting K~ atomic data. 



IV. PRODUCTION OF A NUCLEAR BOUND STATES 

K~ nuclear bound states are expected, generally, to have widths of order 
r ps — 2 ImV r opt (po) ~ 100 MeV. However, for particularly deep states, bound by over 100 
MeV, the dominant two-body pionic decay modes of the KN system get blocked so that 
these states stand a chance of being resolved in some appropriate production reactions. If 
such relatively narrow deeply bound states are ever observed, then the K~ nucleus optical 
potential must be sufficiently deep. Kishimoto [19|] has suggested to search for deeply bound 
K~ nuclear states using the forward (K~,p) reaction in which the incoming K~ meson 
knocks out in the forward direction a bound proton from the target nucleus, while itself get- 
ting captured in a nuclear bound state generated by V opt . This is equivalent to a backward 
(6 = 7r) K~p elastic scattering, and will be denoted by K~p — > pK~ . In this section we 



calculate, more rigorously than done in Ref. |19| , the relevant production cross sections at 
laboratory momentum p-^ = 1 GeV/c. 

The forward differential (K~,p) laboratory cross section on a nuclear target, due to a 
single-particle transition from the proton (n p , l p ) shell to the K~ (n K -, l K -) bound state in 
the residual nucleus, is expressed in the distorted-wave (DW) impulse approximation [|30 
terms of the Fermi-averaged forward K~p — > pK~ laboratory cross section: 



in 




P. 

Iplp >S *' — 1 



/Ax(0 c 

\dn^ 



Here, a is a kinematical factor which accounts for the transformation from the two-body 
laboratory system to the many-body laboratory system |37J and Pn^-^ n _i is an effective 
proton number for the transition n p l p — > uk-Ik- given by 

PnJ^n K ^ = t^tt £ £ % | ((nlm) K - | X <->»x£2(r) | (nlm) p ) jp | 2 . 
[zip -+- i) mK _ mp jp=lp±1/2 

(23) 

In Eq. (|23|) , Sj is the j p proton pickup spectroscopic factor in the target (with a maximum 
value of (2j p + 1)) and the suffix j p attached to the matrix element stands for a possible 
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jp dependence of the bound proton radial wavefunction. We note that Eqs. (|22|j2~3|) are 
equivalent to Eq. (3.15) of Ref. for the (K~, 7r~) reaction. For the distorted waves x( r ) 
we use the eikonal approximation, retaining only the exponential attenuation factor, 



Xy (r) X ( K>(r) « ex V (iqz) exp [-^T(b)j , (24) 

where q is the momentum transfer, purely longitudinal at 0°, with a denoting an average 
KN and pN total cross section, and where the nuclear thickness function T(b) is defined by 

/OO P 
p{r) dz ( / T(b) d 2 b = A) . (25) 
-oo J 

Here b is the impact-parameter coordinate in the plane perpendicular to the direction of 
the forward momentum transfer q. The functions T(b) were evaluated numerically, using 
realistic density distributions for p(r). 

We have calculated the DW expression fl23|) for the forward (K~ , p) reaction at incoming 
momentum p^ — 1 GeV/c on 12 C and 28 Si, to the K~ Is state generated by the K~ nucleus 
DD optical potential. The imaginary part of V op t was reduced to 12% of its nominal strength 
in order to account for the reduced phase space for K~ absorption, assuming no change in the 
in-medium properties of the decay products. The results given below are rather insensitive 
to the precise amount of this reduction. By approximating the K~ and proton bound-state 
wavefunctions by harmonic oscillator (HO) wavef unctions, the matrix elements in Eq. ( |2"3"| ) 
reduce to a one dimensional numerical integration in the variable b, as follows: 

= (spj) s ^ l&y ex p {-\ {iq)2 ) 1 Go{cr) |2 (26) 

for 12 C, and 

1 / 1 ,7 xo\ , „ , , „ , , ,7 ,_s ,2 



JO 



P^u K _ = [^-p) S d5/2 - exp [-^(bq) 2 ) | G 2 (a) - G (a)(l - (bqY/2) \ (27) 

for 28 Si, with the distorted-wave integrals 

G k (a) = J™ 2t k+1 exp(-t 2 ) exp (-^T(bt) \ dt , G k (a = 0) = 1 , k = 0, 2 . (28) 

The mean HO size parameter b is defined by 



11/1 




(29) 



in terms of the K~ and proton HO size parameters. 

The accuracy of this approximation is estimated to incur errors of up to 10% in the 
calculated cross sections. In order to provide a concrete and useful check, we have also 
calculated the forward (ti + ,K + ) reaction cross sections at pi = 1.04 GeV/c on the same 
targets as above, leading to the known Is A hypernuclear ground states in the residual 
nuclei. The momentum transfer characterizing this reaction is quite similar to that in the 
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(K~,p) reaction. The calculation of cross sections for the production of A Is bound states is 
identically the same as for the K~ Is bound states, upon the replacement of bound protons 
in Eqs. ( p2| , |2"3"D by bound neutrons. The results of this calculation are shown in Table [VT|, 
with b\ and b n standing for the A and neutron HO size parameters, respectively, where the 
size parameter b is identified by the exp(— r 2 /2b 2 ) exponential factor of the HO wavefunction. 
The ratio of effective neutron number for the DW calculation, A^ DW , to the effective neutron 
number for the plane- wave (PW) calculation, A^ PW , is also given, as a measure of the effect 
of absorption. The input parameters a (da(0°) / dQ^,) = 0.5 mb/sr for 7i + n — > K + A, a = 27.5 
mb and Sj n = 2j n + 1 for the valent p 3 / 2 and d 5 / 2 neutron orbits in 12 C and 28 Si, respectively, 
are the same as in the eikonal calculation of Motoba et al. |39|| , which did not use the HO 



approximation, for the (7r + , K + ) reaction. The agreement between the latter calculation 
and the present one, as shown in Table [VT|, is very reasonable, inspiring confidence in the 
present results for the (K~ ,p) reaction (Table [VI I|) . We note that the preliminary value for 
the (7r + , K + ) forward cross section on 12 C in experiment E336 at KEK jH]] is 15 ± 1 /ib/sr. 

The forward (K~,p) differential cross sections listed in Table |VI]| were calculated as- 
suming a value of 3.6 mb/sr for the factor a (da(0°)/dQ^) in front of P DW on the r.h.s. 
of Eq. fl2~2"|). This value is based on the measured two-body K~p cm. backward cross 
section 1.7 ± 0.1 mb/sr at p L = 1 GeV/c on the value a = 0.69 appropriate to this 
kinematics, and on an estimated Fermi-average reduction factor of 0.59 due to the peaking 
of the two-body elastic K~p backward cross section at 1 GeV/c as function of the incoming 
momentum, following the procedure outlined in Ref. ||38|| . For the distortion, using the value 



a = 40 mb jyj, the suppressive effect on the forward (K~,p) cross sections is considerably 
stronger than for the (tt + ,K + ) reaction. This is partly due also to the extremely small 
spatial extension of the K~ wavefunction which restricts the (K~,p) reaction to the denser 
nuclear region where absorption prevails. In contrast, the extension of the A wavefunction 
in the K + ) reaction is substantially larger. The effect of absorption gets stronger with 
the atomic number of the target nucleus, as seen clearly from Table |V11| by comparing the 
results for Si with those for C. Another reason for the fast decrease of the calculated cross 
section with increasing A is the gradual increase of the momentum transfer q in the (K~ ,p) 
reaction due to the increased binding of the K~ meson. This is just opposite to the trend 
observed in Table |VI] for the K + ) reaction, where q gradually decreases as function of 
A due to the increased binding of the A hyperon. 

Finally, we comment on the marked disagreement between the results of the present 
calculations and the estimates due to Kishimoto |H| shown in Table |V11[ The wide range 
of estimated cross-section values in Ref. |19[ reflects primarily the dependence of the PW 
expressions Eqs. (^,^) on the HO size parameter b%- ■ The smaller b^- is, the higher the 
PW cross section becomes. However, the lowest value of b K - assumed there, bx- = &at/8 1//4 , 
was erroneously chosen instead of the considerably larger value b K - = &JV/2 1 / 4 that the 
correct HO scaling argument leads to. This should narrow appreciably the range of values 
for the PW cross section. Furthermore, it appears inconceivable that just one value for 
the distortion factor p DW /p PW ; as quoted from Ref. [T!| in Table |V1J| , can be considered 
representative for the whole range of bx~ values assumed by Kishimoto. For this matter, 
the smaller bx- is, the more suppressive is the effect of the distortion. The dependence on 
bx- is rather strong and thus largely cancels the opposite trend of the PW calculation as 
explained above. The overall dependence on bx- is therefore much weaker than one is led to 
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believe in Ref. [0|. As for the absolute scale of the suppression provided by the distortion 
factors p DW /p PW ; we find it inconceivable that Kishimoto's suppression effects are even 
weaker than we and other works |39|] find for the (tt + , K + ) reaction. The present results call 
for extreme caution when contemplating (K~,p) experiments aimed at identifying deeply 
bound K~ nuclear states. 



V. SUMMARY AND CONCLUSIONS 

At present the best tool for exploring the K~ interactions in the nuclear medium at 
low energy is the study of strong interaction effects in kaonic atoms. Extrapolating this 
interaction to higher densities, as encountered in astrophysical scenarios, must rely on some 
theory, beyond the phenomenological potential that fits the data very well. The prime aim 
of the present work was to see how far one can go in bringing microscopic approaches to the 
A^-nucleus interaction into agreement with the K~ atomic data. The chirally motivated 
coupled channel approach to the KN interaction, which is quite successful in reproducing 
all the low energy K~p data, was chosen as a starting point. Earlier attempts to use 
this approach indicated poor agreement with the atomic data and only by introducing self 
consistency into the theory it became possible to achieve barely acceptable fits. Empirical 
modifications of the interaction managed to improve the fits to the atomic data, but at the 
cost of losing contact with the underlying KN interaction model, thus making questionable 
any extrapolation of such optical potentials to higher densities. 

In the present work we addressed the problem by requiring simultaneous fits to atomic 
and K~p data within the chirally motivated coupled channel approach, including self con- 
sistently the self energies of the kaons and, in some cases, of the nucleons. It was found that 
minor modifications of parameters of this theory led to good agreement with the atomic 
data while maintaining the good agreement with the K~p data. In fact, we found that one 
can fit only the atomic data and still not lose the good fit to the K~p data. The best value 
of x 2 P er point obtained here for the K~ atomic data is 2.2, compared to 1.5 for the best-fit 
phenomenological potential. Nevertheless, the fits of the self consistent chirally motivated 
coupled channel potentials are quite good. 

The depth of the real part of the K~ nucleus potential is also of prime interest in 
connection with the possibility of kaon condensation in collapsing stars. Deep (180 MeV) 
potentials were found in the phenomenological analysis, whereas shallow (55 MeV) potentials 
are found in the framework of models which require the K~ optical potential to be derived 
self consistently. It seems impossible to reconcile the phenomenological potentials with the 
present microscopic potentials both in terms of depth and in terms of values of \ 2 - Since 
the depth of the A'~-nucleus optical potential cannot be resolved by studying only kaonic 
atoms, we briefly discussed the ability to do so with the help of reactions initiated by stopped 
K~ mesons. The mechanism behind it is the sensitivity of the DW integrals to the depth 
of the real part of the K~ potential, due to the node structure of the wavefunction, which 
is quite different from each other for deep and for shallow potentials. It was demonstrated 
that A-hypernuclear ground-state production rates calculated for the {K^ op ,n) reactions on 
carbon differ by more than a factor of 3 between the different potential depths mentioned 
above. Unfortunately, these calculated rates are still several times smaller than the measured 
ones P7[ , so that these reactions cannot yet be used to reach a definite conclusion on the 
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depth of the K~ nucleus optical potential near threshold. However, it is plausible that by 
studying the effect of the K~ potential on the pion spectrum of the capture reactions, one 
could reach more definite conclusions. 

Finally, we addressed the possibility of identifying deeply bound K~ nuclear states in 
the forward (K~,p) reaction on nuclear targets. Using the very deep DD potential ||, we 
calculated a cross section as large as 47 /xb/sr at p-^ = 1 GeV/c for the production of the 
Is K~ bound state on 12 C. This K~ nuclear 'ground state' is bound by about 120 MeV, 
which strips it off most of the phase space for its dominant pionic decay modes, so that 
its residual width could become as small as 10-20 MeV. In contrast, using the HO size 
parameter h K - = 1.22 fm from Table |V11| , the first excited (lp) state is estimated to lie 
about hu ~ 50 MeV higher and, therefore, it should be much wider and hardly observable. 
We verified this estimate by an explicit calculation. If the K~ nuclear potential is not of the 
'deep' kind, it may still accommodate bound states, but these will be very wide, judging from 
the depth of Im V opt found in this work for the 'shallow' chirally motivated potentials. We 
calculated the ls K - bound state production cross section also on 28 Si, finding it to be about 
8 times smaller than on 12 C. This rules out using medium- weight or heavy targets for this 
reaction. Not much different production cross sections should be expected at lower incoming 
momentum, say at 600 MeV/c. Our calculated cross sections, tested for the kinematically 
similar (n + ,K + ) reaction, are substantially lower than those estimated by Kishimoto [|19|j. 
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TABLES 



TABLE I. Kaonic atom optical potentials using two sets of self consistent in-medium KN 
scattering amplitudes, with or without adding a complex l tp' potential, and with or without scaling 
the 1 = 1 amplitude. Vr, and V\ refer to the depth of V op t for Ni. 



amplitude 


t R (fm) 


ti(fin) 


sr 


si 


X 2 /N 


T/ R (MeV) 


Vi(MeV) 


Ref. | 


15 


i 








1 


1 


4.62 


-44 


-54 








0.11±0.02 


-0.20±0.02 


1 


1 


2.73 


-58 


-23 














2.0±0.2 


0.9±0.2 


3.08 


-61 


-52 








1.17±0.11 


-0.03±0.15 


-0.5±0.5 


0.96±0.54 


2.00 


-170 


-49 


Ref. | 


m 


i 








1 


1 


12.7 


-34 


-62 








1.0±0.1 


0.26±0.06 


1 


1 


2.46 


-159 


-95 














4.5±0.3 


5.7±0.7 


2.82 


-134 


-116 








2.9±0.4 


-2.8±0.8 


-6.1±1.7 


34±9 


1.49 


-204 


-82 



TABLE II. Calculated branching ratios fll3| ) at the K p threshold, and total cross sections (in 
mb) for selected K~p reactions ( |i~4l) at incident momentum pl = HO MeV/c, for the 'no SC and 



Vn, full' parameterizations of Table III. Also shown are the experimental data. 







7 




R n 


a(K p) 


a(K°n) 


ct(7t+£-) 


o-(tt-£ + ) 


'no 


SC 


2.383 


0.667 


0.171 


93.2 


41.2 


60.1 


32.3 


'V N , 


full' 


2.347 


0.669 


0.198 


94.5 


38.0 


71.2 


42.0 



exp. 



2.36±0.04 0.664±0.011 0.189±0.015 92 ± 8 29 ± 6 64 ± 10 29 ± 6 



TABLE III. Values of x 2 > optical potentials and scaling factors for the present chiral model. 
Xatom refer to the K~ atom 65 data points and \\ is the sum of the former and the x 2 for the 
7 data points for the free-space K~p data. Vrj are the optical potential values at the center of 
the Ni nucleus and Vn stands for fits where the nucleon optical potential is included. The scaling 
factors / and ct(j) multiply the parameters f n and ctj, respectively, of the coupled channel model. 



model 



Xt Xltom Vk(MeV) Va(MeV) / 



a(K-p) a(K°n) a(vr°A) a(7rS) 



no SC 
SC 

atoms 
full 



1082 1076 

443.1 436.7 

169.9 157.2 

165.3 157.5 



Vn, atoms 156.3 145.4 
V N , full 154.0 145.8 



-117 -67 11111 

-70.6 -85.6 1 1 1 1 1 

-66.9 -55.2 1.12±0.05 1.00L0.08 1.21±0.08 1.05±0.17 1.39±0.07 

-66.8 -54.0 1.14±0.05 1.03±0.07 1.24±0.08 1.13±0.16 1.35±0.07 

-56.3 -61.8 1.04±0.04 0.86±0.06 1.07±0.05 1.00±0.14 1.18±0.08 

-54.9 -60.2 1.06±0.04 0.90±0.06 1.10±0.06 1.08±0.15 1.17±0.07 
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TABLE IV. K optical potentials and calculated strong-interaction shifts and widths, in eV, 



for the 1p and 3d K atomic levels in C. 



Vopt 


a (fm) 


B (fm) 


7 


£2p 


^2p 


e3d 


^3d 


chiral 








-616 


1484 


+0.01 


0.57 


t e s 


0.63 + i0.89 






-590 


1360 


-0.01 


0.64 


teS 


1.30 + i0.80 






-605 


1731 


-0.01 


0.85 


DD 


-0.15 + i0.62 


1.65 - i0.06 


0.23 


-468 


1578 


-0.08 


0.64 


exp. 35[ 








-590 ± 80 


1730 ± 150 




0.98 ± 19 



TABLE V. Calculated capture rates on 12 C per stopped K (in units of 10 3 ) to the summed 
Pn — ► s\ 1~ excitations in j^C and ^ 2 B, for the K~ optical potentials of Table |TV| . 

final chiral t e g t e fi DD 

fC 0.231 0.169 0.089 0.063 

1 2 B 0.119 0.087 0.046 0.032 



TABLE VI. Input and results of DWIA eikonal calculations for the forward (-7r + , K + ) reaction 
cross section (in /ib/sr) at incoming momentum pl = 1-04 GeV/c, see text for details. 



target 


R A 


q 


b A 


K 


N BW /N PW 


(d ( r(0°)/ ( ifi L )^ + ^ + ) 


nucleus 


(MeV) 


(MeV/c) 


(fm) 


(fm) 


[present] 


[present] Ref. || 


12 C 


10.8 


335 


1.72 


1.52 


0.237 


15.4 17.4 


28 Si 


16.7 


322 


1.92 


1.72 


0.132 


7.7 8.9 



TABLE VII. Input and results of DWIA eikonal calculations for the forward {K , p) reaction 
cross section (in ^ib/sr) at incoming momentum = 1 GeV/c, see text for details. 



target 


B ls 


Q 


b K - 


Go 


G 2 


pDW /pPW 




nucleus 


(MeV) 


(MeV/c) 


(fm) 






[present] Ref. || 


[present] Ref. @ 


12 C 


122 


369 


1.22 


0.308 




0.095 0.25 


47 100-490 


28 Si 


144 


404 


1.42 


0.183 


0.269 


0.040 0.16 


6.0 35-180 
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FIG. 1. Cross sections for K~p scattering and reactions to the channels indicated in the figure. 
Results of the full fit ('Vn, full', see Table III), displayed by solid lines, are compared with the 
available data (see text). Results for the free-space chiral-model parameterization ('no SC) are 
shown for comparison (dashed lines). 
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FIG. 2. Real (top) and imaginary (bottom) parts of the isospin-averaged K~N (effective) 
threshold scattering amplitude as function of density p/po, calculated for the 'no SC (dashed 
line), 'SC (dot-dashed line) and 'Vn, full' (solid line) chiral-model parameterizations (see Table 
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FIG. 3. The KN free-space scattering amplitude fi=o as function of cm. energy. Results are 
presented for the 'no SC (dashed lines) and 'Vn, full' (solid lines) chiral-model parameterizations 



from Table III 
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